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Abstract 

We present preliminary results of searches for exclusive charged-i? decays to K^tt^tt^ from 
61.6 million BB pairs collected at the T(4S) resonance with the BABAR detector at the SLAC PEP- 
II asymmetric B Factory. The Dalitz plot is divided into eight regions and, using a maximum- 
likelihood fit, we measure statistically significant yields in all regions. We interpret the results as the 
following branching fractions averaged over charged-conjugate states: B{B^ K*^ (S92)t^^ , K*^ — >■ 
K+TT-) = (10.3 ± 1.2l^:?) X 10-^ B{B+ fo{980)K+,fo tt+tt") = (9.2 ± 1.2l|^) x 10-^ 
B{B+ XcoK+, Xco TT+TT-) = (1.46 ± 0.35 ± 0.12) x 10"^ and B{B+ D°tt+, D° K+tt-) = 
(184.6 lb 3.2 lb 9.7) x 10~®. The first uncertainty is statistical and the second is systematic and 
includes resonance-model and interference uncertainties. We give 90% confidence- level upper lim- 
its on the branching fractions of the following channels: B{B~^ p^{770)K~^) < 6.2 x 10~^ and 
B{B+ K+TT-TT+ non-resonant) < 17 x 10"*^. 
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1 Introduction 



The study of charmless hadronic B decays can make important contributions to the under- 
standing of CP violation in the Standard Model as well as to models of hadronic decays. The 
measurement of ^ decays to the final state K'^ir'-ir'^ via intermediate resonances can be used to 
search for direct CP violation. The three-body final state is unique in the search for weak phases 
since it is possible to isolate the strong phase variation for overlapping resonances. There has been 
recent theoretical progress on proposed methods for extracting the Cabibbo-Kobayashi-Maskawa 
angle 7 through the interference of XcoK~^ with other K^tt~tt~^ final states 0l2l- Study of these 
decays can also help clarify the nature of the resonances involved, not all of which are well under- 
stood. We present preliminary results on the branching fractions of decays to the final state 
K~^TT~TT~^ both non-resonant and by way of intermediate resonances. 

2 The BaBAR Detector and Dataset 

The data used in this analysis were collected with the BABAR detector at the PEP-II asymmetric- 
energy e~^e~ storage ring at SLAC. The data sample consists of 61.6 million BB pairs, correspond- 
ing to an integrated luminosity of 56.4 fb~^ collected on the T(4S) resonance (10.58 GeV) during 
the 2000-2001 run. In addition, a total integrated luminosity of 6.4 fb~^ was taken 40 MeV below 
the T(4S) resonance, and was used to study backgrounds from continuum production. 

The BABAR detector is described in detail elsewhere 0; the main parts relevant for this analysis 
are the tracking and particle identification sub-systems. 

The 5-layer double-sided silicon vertex tracker (SVT) measures the impact parameters, angles, 
and transverse momenta of tracks. Outside the SVT is a 40-layer drift chamber (DCH), which 
measures the transverse momenta of tracks from their curvature in the 1.5-T solenoidal magnetic 
field. The ionization energy loss of charged tracks, dE/dx, in the SVT and DCH is used in the 
particle-type identification. The tracking system has a momentum resolution of 0.5% for a trans- 
verse momentum of 1.0 GeV/c and a typical dE/dx resolution of 7.5%. 

Surrounding the DCH is a detector of internally-reflected Cherenkov radiation (DIRC), which 
provides charged-hadron identification in the barrel region. The separation between pions and 
kaons varies from > 8a" at 2.0GeV/c to 2.5a at 4.0GeV/c, where a is the average resolution on the 
Cherenkov angle. 

The DIRC is surrounded by a Cesium Iodide electromagnetic calorimeter (EMC), which is used 
to measure the energies and angular positions of photons and electrons with excellent resolution. 
In this analysis, the EMC is used to veto electrons. 

3 Analysis Method 

The K^TT^TT'^ final state can be represented in a Dalitz plot 1^. The many resonant B decay 
modes form bands in such a plot. These resonances often overlap and interfere so the whole Dalitz 
plot should be considered before assigning a branching fraction to a specific mode. This analysis 
divides the Dalitz plot into regions, each of which is expected to be dominated by a particular 
contribution. First the yields in these regions are determined, using a maximum-likelihood fit, 
with no assumption on the form of the intermediate resonances. We then interpret these yields as 

^Charge-conjugate states are included throughout this document. 
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branching fractions, assuming a model for the contributions to the Dahtz pfot. We also consider 
the uncertainty of the model and the effect of overlap and interference between these contributions. 

3.1 Dalitz Plot Regions 

The K^TT^TT^ Dalitz plot is divided into eight regions. Each region is designed to contain a large 
proportion of the decays of the expected dominant resonance (if any) and to minimize contributions 
from neighboring modes. The definitions of the regions are given in Table ^^-nd they are illustrated 
in Figure n 

Regions I, II and III are characterized by narrow bands in the invariant mass of the K~^tt~ 
system, rnxn- Region I is expected to be dominated by the K*^{892). The primary resonance 
contribution to region II, labeled "higher K*^" , is not currently known. The areas where these 
bands cross the vrvr resonances {p{770), /o(980), "higher /" and Xco) are excluded to limit biases from 
interference. Region III is dominated by the production of D^tt^. The relatively high branching 
fraction for this mode allows it to be used to correct for differences between data and simulated 
events and evaluate systematic uncertainties. 

Regions IV, V and VI are characterized by narrow bands in the vr+vr^ invariant mass, m^Tf. 
Regions IV and V are expected to be dominated by the p{770) and /o(980), respectively. The 
resonance contributions to region VI ("higher /") are not well defined. The area where these 
regions would intersect the band, 1.8 < rriKw < 1.9GeV/c'^, is excluded from regions IV, V and 
VI. The area where the other Kir resonances cross is not excluded from regions IV, V and VI as the 
overall interference uncertainty on B{B^ — > /o(980)-Rr^) and B{B^ p^{770)K^) was estimated 
to be smaller when this area is excluded. Region VII, denoted "high mass", could contain higher 
charmless and charmonium resonances as well as a non-resonant contribution. 

Region VIII is dominated by XcoK^. A lower limit on rriKn ensures that this region is free of 
contamination from resonances in regions I, II and III. The area 3.355 < rriT^T^ < 3.475 GeV/c^ is 
removed from all other regions to avoid this charmonium background. 

Table 1: Regions of the K~^TT~7r~^ Dalitz plot. The regions are kinematically defined by the m^Tr 
and m-TTTr mass selection criteria. The symbol "!" is used to imply exclusion. The assumed dominant 
contribution is noted for each region. 



Region 


Dominant 
Contribution 


Selection Criteria 


TJiKn (GeV/c^) 


m^TT (GeV/c^) 


I 


K*"(892)7r+ 


0.816 < niK^ < 0.976 


m^^ > 1.5 !(3.355 < m^^ < 3.475) 


II 


"higher 7r+ 


0.976 < ruK^ < 1.8 


m^^ > 1.5 !(3.355 < m^^ < 3.475) 


III 




1.835 < rriK^ < 1.895 


!(3.355 < m^^ < 3.475) 


IV 


p^{770)K+ 


!(1.8 < ruKn < 1-9) 


0.6 < rn-TT^ < 0.9 


V 


fo{980)K+ 


!(1.8 < rriKn < 1-9) 


0.9 < rn-TT^ < 1.1 


VI 


"higher /" K+ 


!(1.8 < ruK^ < 1.9) 


1.1 < TTlTTTr < 1.5 


VII 


"high mass" 


mK-K > 1-9 


m^^ > 1.5 !(3.355 < m^^ < 3.475) 


VIII 




mK-K > 1-9 


3.37 < m^^ < 3.46 
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Figure 1: An illustration of the different regions in the Dalitz plot as defined in Tabled 



3.2 Candidate Selection 

B candidates are reconstructed from charged tracks that have at least 12 hits in the DCH, a 
maximum momentum of 10 GeV/c, a minimum transverse momentum of 100 MeV/c, and originate 
from the beam-spot. The B candidates are formed from three-charged-track combinations and 
particle identification selections are applied. Mass hypotheses are assigned accordingly and the B 
candidates' energies and momenta are required to satisfy appropriate kinematic constraints. 

Two kinematic variables are defined that are included in the maxim um-likeliho od fit described 
later. The first of these is the beam-energy-substituted mass tues = \l {E^ — The energy of 

the B candidate is defined as E}, = (^s + po •P_b)/-E'o, where ^/s and Eq are the total energies of the 
e^e~ system in the center-of-mass (CM) and laboratory frames, respectively, and po and p^ are 
the three-momenta in the laboratory frame of the e^e~ system and the B candidate, respectively. 
The mES value should be close to the nominal B mass for signal events. 

The second variable used is the energy difference, Ai?, between the energy in the CM of the 
reconstructed B candidate, E*^, and the beam energy, /S.E = E*^ — \/i/2. /S.E is dependent on 
the mass hypotheses of the tracks. To each track a mass is assigned appropriate for the particle 
identification selections applied. For signal events, /^E should be centered at zero. 

To identify charged pions and kaons, we use dE/dx information from the SVT and DCH 
for tracks with momenta below 700 MeV/c, the number of photons measured by the DIRC for 
tracks with momenta above 500 MeV/c, and the Cherenkov angle for tracks with momenta above 
700 MeV/c. Kaons are selected with requirements on the product of likelihood ratios determined 
from these measurements and pions are required to fail the kaon selection. The average selection 
efficiency for kaons in our final state that have passed the tracking requirements is ~ 80% including 
geometrical acceptance, while the misidentification probability of pions as kaons is below 5% at all 
momenta. The kaon veto on pions in our final state is ~ 98% efficient. We veto electron candidates 
by requiring that they fail a selection based on information from dE/dx, shower shapes in the 
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EMC and the ratio of the shower energy and track momentum. The probabihty of misidentifying 
electrons as pions is approximately 5%, while the probability of misidentifying pions as electrons is 
- 0.2%. 

3.3 Background Suppression and Characterization 

The dominant background in this analysis is from light quark and charm continuum production. 
This background is suppressed by imposing requirements on topological event shape variables, and 
the remainder is characterized and parameterized in the maximum-likelihood fits used to extract the 
signal yield from the data. There are also backgrounds from B decays which, although contributing 
far fewer events, can be more difficult to separate from the signal and must also be parameterized 
in the fit, which is described in detail later. 

The event shape variables used to suppress continuum background are calculated in the T(4S) 
rest frame. The first is the cosine of the angle 9t between the thrust axis of the selected B candidate 
and the thrust axis of the rest of the event, i.e., of all charged tracks and neutral particles not in 
the B meson candidate. For continuum backgrounds, the directions of the two axes tend to be 
aligned because the daughters of the reconstructed candidate generally lie along the dijet axis of 
such events so the distribution of cos 9t is strongly peaked towards ±1. For B events, however, 
the distribution is isotropic because the decay products from the two B mesons are independent 
of each other and the B mesons have very low momenta in the T(4S) rest frame. To improve the 
signal-to-background ratio, the criterion |cos0t| < 0.9 is applied to all regions except region III 
(D^). This selection removes about 60% of the continuum background while retaining over 90% of 
the signal. The cos 9t requirement in the region is varied to estimate a systematic uncertainty 
on the efficiency of this criterion in the other regions. 

We also make use of a Fisher discriminant ,5i, using a linear combination of the angle between 
the B candidate momentum and the beam direction; the angle between the B candidate thrust 
axis and the beam direction; and the energy flow of the rest of the event into each of 9 independent 
concentric 10° cones around the thrust axis of the reconstructed B [S]. The variables and weights of 
the Fisher discriminant were chosen to optimize the separation of our final state and the continuum 
background after the cos 9t criterion has been applied. The resulting Fisher variable, is used in 
the maximum-likelihood fit. 

The B decay backgrounds are from four main sources: combinatorial background from three 
unrelated tracks; three- and four-body B ^ D decays; charmless four-body decays with a missing 
track and three-body decays with one or more particles misidentified. These backgrounds are 
reduced by the particle identification selections and, where possible, removed by vetoing regions of 
the invariant-mass spectra of pairs of the final-state particles. The influence of remaining specific 
backgrounds on the signal yield obtained from the maximum-likelihood fit was established using 
test fits with Monte-Carlo simulated data (MC) with the expected number of signal, continuum 
background and B background events. Background modes that significantly contributed to the 
signal yields in these tests are parameterized for the final fit to the data. However, if the background 
contributes only a few events, it is instead subtracted from the signal yield. 

The combinatorial background from B decays is less than 2% of the continuum background. 
The shape of its A.E and mss distributions are similar to the continuum background and these 
events were found to be fitted as such in the test fits so no additional parameterisation was required. 

The particle-misidentification background has several sources. B^ — > J/tpK^ and B^ — > 
^{2S)K~^ decays which contribute through muon/pion misidentification are removed completely 
by excluding events with 2.97 < 771^^^^ < 3.17GeV/c^ or 3.56 < m-n-T^ < 3.76GeV/c^. This also 
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excludes other B backgrounds containing J/il^ or ^{2S), such as B'^ — > J/ij^K*'^ . The contributions 
from B^ — > K^e^e~ are neghgible due to the electron veto. Non-resonant B~^ — > K~^ii~^n~ decays 
contribute with 18ib7 events uniformly distributed over the Dalitz plot, using the branching fraction 
from Events are expected only in regions II, VI and VII where they are subtracted from the 
yields. Contributions from the B^ — > 7r"'"7r~7r"'" and B^ K^K~tt~^ final states may affect regions 
I, II, IV, VI and VII. The channels concerned have not yet been observed and so an additional 
negative uncertainty is added to the signal yield equal to the number of events expected in the fit 
corresponding to the upper limit measured in JH] . This is 20 events for region VII and 6 events or 
fewer for all other regions. 

The decay B~^ — > D^tt^ A'+vr^vr^ contributes to the measured signal yield when not 
parameterized in the test fits in regions II and VII, while B^ (770) with K^'k~ 

and p"*" -K^TT^ has a significant effect on the test fits in region VII. These modes are therefore 
parameterized in the final fit for the affected regions. 

The decay B~^ — > r]'K~^ with rj' — > p^{770)"f, vr^vr^ is the only charmless channel with a 

four-body final state found to contribute significantly. The expected number of events is 31 and 12 
in regions IV and V, respectively, using the branching ratio measured in [S]. If not parameterized, 
these contribute to the signal yields in the test fits. This mode is therefore included in the final fit 
for those regions. 

3.4 The Maximum-Likelihood Fit 

We form Probability Density Functions (PDFs) for the 3 variables ttt-es; and in each 
region. For each hypothesis / (signal, continuum background, and if applicable, B background) 
these three PDFs form a product Pi, which models that particular hypothesis. These products 
are functions of the variables x and parameters a of the PDFs. The likelihood for an event j 
is formed by summing the products over the M hypotheses, with each product weighted by the 
number of events in that hypothesis n;. A product over the N events in the data sample of the 
per-event likelihoods along with a Poisson factor forms the total likelihood function £, written in 
equation (Q). 

(M \ N / M \ 
i=i / j=i \i=i J 

The fit is performed in two stages for each region. First, one-dimensional fits are performed 
on the particular data samples detailed below in order to determine the PDF parameters. Then 
the multivariate fit is performed on the final data samples to extract the signal and continuum 
background yields. 

The signal PDF parameters, particularly the width of the Ai? distribution, vary across the 
Dalitz plane. Therefore, they are found for each region separately using Monte-Carlo-simulated 
signal of the expected dominant resonance where available and, otherwise, non-resonant K~^ti~ii~^ 
selected for that region. Some differences have been observed between MC and data in the dis- 
tributions of the mESi and T variables. These differences are measured in the high-statistics 
B^ D^TT^ dominated region and are used to correct all regions where all the signal PDF param- 
eters are fixed in the final fits. The mss PDF is a Crystal Ball function ^Uj, the AE PDF is two 
Gaussians with equal means and the PDF is a "bifurcated" Gaussian (a Gaussian with different 
widths above and below the mean). 
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The background PDF parameters are also found to vary across the Dalitz plane and so are 
determined individually for each region. The T PDF is a bifurcated Gaussian where the parameters 
are determined using data sidebands defined by 0.1 < |Ai?| < 0.35 GeV except in regions II and 
VII where, due to B background, only the positive A£' sideband is used. The rwES variable 
is parameterized by the Argus threshold function |llj with two parameters: a fixed kinematic 
endpoint and a shape parameter which is left to float in the final fits. The PDF is a first-order 
polynomial where the gradient is also left to float in the final fits. 

4 Physics Results 
4.1 Fit Results 

The signal yields for the various regions of the /C+7r^7r+ Dalitz plot are shown in Table |21 The 
first uncertainty is statistical and the second is systematic. The systematic uncertainty is from 
the uncertainties on the PDF parameters and from the B background subtraction. The yields are 
found to be statistically significant in all regions (> 5cr), where the statistical significance is taken 

as ^-2 log (£max / {n^^g =0) ) • 

Table 2: Signal yields for the regions of the Dalitz plot before and after B background subtraction. 
The first uncertainty is statistical, the second systematic. 



Region 


Signal Yield 


Signal Yield after B 
background subtraction 


I 

II 

III 

IV 

V 

VI 

VII 

VIII 


161 ±18 ±4 
406 ± 28 ± 12 
3755 ± 66 ± 11 
66 ± 15 ± 3 
179 ± 19 ± 5 
127 ± 19 ± 5 
147 ± 23 ± 7 
26 ± 6 ± 1 


161 ±18 ±4 
405 ± 28t}^ 
3755 ±66 ±11 
66 ± 15+? 
179 ±19 ±5 
126 ±19 ±5 

133 ± 23t^2 
26 ± 6 ± 1 



The projection plots of rriES) Ai? and T for each region are shown in Figures [21 and IHI To 
produce these plots the projected variable is excluded from the likelihood functions and the ratio of 
the signal and background likelihoods for each event evaluated. The histograms are produced with 
a selection on this per-event likelihood ratio, where the selection value was chosen separately for 
each region to best illustrate the signal contribution. The projection of the fits onto that variable 
is then superimposed. 

Figure shows the Dalitz plot for data events within the signal region 5.2715 < txies < 
5.2865 GeV/c^ that have a per-event likelihood ratio, formed from the AE and PDFs, greater 
than 5. To illustrate the expected background distribution, events passing the same likelihood 
selection but having a value of mES between 5.25 < tues < 5.26 GeV/ are shown alongside. The 
size of this sideband is chosen to contain approximately the same number of background events as 
are in the signal region. The mass intervals close to J/^ and ^{2S) are removed. 
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Figure 2: Projection plots in mESj ami produced by selecting on the event likelihood 
ratio formed from the other two fit variables for (from top to bottom) regions I, II, III and IV. 
The superimposed curve is a projection of the full fit with the background component shown as a 
dotted line and, for the signal component shown as a dashed line. 
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Figure 3: Projection plots in rriESi and produced by selecting on the event likelihood ratio 
formed from the other two fit variables for (from top to bottom) regions V, VI, VII and VIII. The 
superimposed curve is a projection of the full fit with the background component shown as a dotted 
line and, for the signal component shown as a dashed line. 
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Figure 4: Dalitz plots showing the observed distribution in the signal tties region (left) and contin- 
uum background from the mES sideband (right) with the mass intervals close to J/^ and ^{2S) 
removed. 
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Figure 5: Background-subtracted and efficiency-corrected projections of the Dalitz plot in rriKn 
and m-jriT in the ranges 0.6GeV/c^ < rnKw < 1.8GeV/c^ and 0.2GeV/c^ < m^^ < 1.5GeV/c^ with 
J/^, ^{2S) and vetoes apphed. Peaks at the K*^{892) and /o(980) masses are clearly visible. 



Figure IS] shows background-subtracted projections of the Dalitz plot in rriKn from 0.6GeV/c^ to 
1.8 GeV/ and m^rTr from 0.2 GeV/ to 1.5 GeV/ c^. For these plots, signal events are obtained using 
a likelihood ratio selection on the tues find J- PDFs and the requirement |A£^| < 0.07 GeV, while 
the background is determined from the sideband 0.1 < AE < 0.35 GeV. The J/^ and ^(25) 
vetoes are applied. The plots are efficiency-corrected using non-resonant K~^TT~Tr~^ Monte-Carlo 
data. Peaks at the K*°(892) and /o(980) masses are clearly visible. It is not clear what other 
channels are contributing, although there is a large signal in the region 1.1 < mx-K < 1-4 GeV/ c^. 
Figure El shows background-subtracted plots of the resonant mass and helicity angle distributions 
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Figure 6: Projection plots of the resonant mass and cosine of the helicity angle for regions I and V. 
The background has been subtracted from all plots and, for the helicity distributions, a correction 
applied to each bin to account for the variation in reconstruction efficiency. The bold lines are the 
mass and helicity angle distributions that would be expected from the dominant resonant mode, 
i^*''(892)7r"'" for region I and /o(980)i^"'" for region V with Breit-Wigner lineshapes 

and the masses and widths given in Tabled 



for regions I and V, where the 7^**^(892) and /o(980) resonances are expected to dominate. Only the 
helicity distributions have been efficiency-corrected as the efficiency does not vary over the range of 
the plotted mass distributions. The plots in Figure El have been overlaid with the distribution of the 
expected dominant resonance using Breit-Wigner line-shapes for the mass distributions, cos^ 6h for 
the K*^{892) angular distribution and a flat line for the scalar /o(980) angular distribution. There 
is good agreement between the overlaid and observed distributions indicating that the expected 
resonances are indeed dominant in these regions. 

Figure [7| shows plots for region II produced in the same way. It is not clear from these distribu- 
tions which resonances are present but it is unlikely that a single resonance can describe the uik-k 
distribution. Superimposed is the effective-range parameterization used for the scalar i^Q'^(1430) 
observed in the LASS experiment |12|, using parameters taken from ^SI- The agreement between 
this parameterisation and our observed distribution is good up to 1.6 GeV/cP. The increase in event 
density above 1.6 GeV/c^ is most likely caused by contributions from other higher resonances. The 
possibility that the is contributing to this plot has been excluded by MC studies. 
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Figure 7: Projection plots of the resonant mass and cosine of the hehcity angle for region II. The 
background has been subtracted from all plots and, for the helicity distributions, an efficiency 
correction applied to each bin. The mass parameterization obtained from the LASS experiment 
is superimposed as is a flat scalar distribution for the cosine of the helicity angle. The increase 
in event density above 1.6 GeV/c^ is most likely caused by contributions from other higher- mass 
resonances. 



4.2 Branching Fractions 

Regions III and VIII have negligible contributions from K^tt^tt^ channels other than 

D^TT^ and XcqK'^, respectively. D^tt^ and XcqK^ do not contribute to other regions. The branching 
fractions for these channels are defined by: 

where Y is the signal yield and = (61.6 it 0.7) x 10^ is the number of BB events in the sample. 
It is assumed the T(4S) decays equally to neutral and charged B meson pairs. The reconstruction 
efficiency, e, is calculated using signal MC and is corrected for MC/data discrepancies in tracking 
and particle identification. We achieve efficiencies of 0.330 it 0.017 and 0.288 it 0.017 for D^tt^ and 
XcoK~^, respectively, where the uncertainty is purely systematic. 

For the other regions, we calculate the branching fractions from the measured yields taking into 
account the resonance cross-feed and model dependence. Y becomes a vector of the yields in each 
Dalitz region, B becomes a vector of the branching fractions and the efficiency becomes a matrix 
M the elements of which are the probability of an event of a particular decay to be found in a 
particular region. 

B = M-i Y/N^^ (3) 

The branching fractions measured depend on the model of resonances assumed in calculating the 
efficiency matrix M. We split M into two component matrices, P and e, such that each element 
Mjj = Pijeij. The P matrix contains the event distribution around the Dalitz plot and the e matrix 
contains the reconstruction efficiencies, so the model dependence is contained in the P matrix. 

We assume one dominant contribution per region and the decay modes in the chosen resonance 
model are given in Table 01 The i^T*" (892)7r+ and /o(980)i^+ channels have been seen and there is 
evidence for ( 770) jl4l [W\ . For regions II, VI, VII, there are a number of possible contributions: 
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Table 3: A summary of the model used to calcu 



i 


Decay Mode 


Lineshape 


Mass (MeV/c^) 


Width (MeV/c^) 


secondary BP (%) 


Alternative Resonance 


1 
2 
4 
5 
6 
7 


/i:*"(892)7r+ 

/fo*°(1430)7r+ 

p°{770)K+ 

M980)K+ 

f2{l270)K+ 

non-resonant 


BW 

BW (LASSHa) 

Blatt-W 
BW (FlatteUHl) 
BW 
flat 


896.10 ± 0.27 
1412 ± 6 
769.0 ± 0.9 

980 ± 10 
1275 ± 12 

all masses 


50.7 ± 0.6 
294 ± 23 

150.9 ± 1.7 
70 ±30 
185 ± 30 


33.3 
93± 10 
100 
dominant 
85 ±2 


ii:2*°(1430),ii:f' (1680) 
/o(1370), /2(1430) 



ate branching fractions. 



Table 4: The elements of the matrix P used to calculate the branching fraction central values, 
where Pjj is the probability for an event of decay mode i to be produced in region j. 



Decay Mode i 


Region j 


I 


II 


IV 


V 


VI 


VII 


1 K*"(892)7r+ 


0.659 


0.099 


0.000 


0.000 


0.144 


0.000 


2 K*°(1430)7r+ 


0.013 


0.720 


0.019 


0.018 


0.047 


0.110 


4 pO(770)K+ 


0.000 


0.000 


0.738 


0.120 


0.052 


0.000 


5 /o(980)K+ 


0.000 


0.000 


0.084 


0.779 


0.084 


0.000 


6 /2(1270)K+ 


0.008 


0.031 


0.020 


0.064 


0.741 


0.111 


7 high mass 


0.013 


0.127 


0.032 


0.029 


0.074 


0.642 



for our model, we choose Kq^{1A30), /2(1270) and a flat non-resonant K~^7r~Tr~^. The masses and 
widths are taken from the Review of Particle Physics pLS^ , and non-relativistic Breit-Wigner line- 
shapes are used for all channels except for the broad p(770) resonance, where we use a relativistic 
Breit-Wigner lineshape with Blatt-Weisskopf damping [T71. The matrix element Pij gives the 
probability of an event of decay mode i to be produced in region j calculated using this model 
and including angular distributions and phase space. The elements are shown in Table HI There 
are large uncertainties in this model: the dominant resonance is unknown in some regions, and 
there are uncertainties on the masses and widths of the resonances, as well as the choice of line- 
shapes. Alternative resonances, line-shapes and the uncertainties on resonance parameters are 
listed in TableOl These cause uncertainties in P and therefore on the branching fractions. "Model" 
uncertainties on the branching fractions are evaluated that take into consideration all of these 
uncertainties in the model. There are also uncertainties in P due to the interference between the 
resonances and these "interference" uncertainties are also evaluated. 

The element Cij is defined as the number of decay-mode i MC events with a candidate in region j 
passing all the selection criteria divided by the number of decay mode i events generated in region j. 
The e matrix is calculated using resonant signal MC for K*^{892)Tr~^ , i?+ p^{770)K~^ and 

fo(980)K^ , and non-resonant MC for other decay modes. These efficiencies have negligible 
dependence on the resonance line-shape; however, they are affected by the angular distribution. 
Corrections for MC/data discrepancies in particle identification and tracking are applied. The 
efficiencies are shown in Table |5J 

Table El gives the branching fractions produced from the vector of yields with statistical, 
reconstruction-systematic, resonance-model and interference uncertainties. 
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Table 5: The elements of the efficiency matrix e, where Cij is the efficiency of decay mode i events 
produced in region j to pass the selection cuts. The statistical uncertainty on these numbers is 
negligible and the fractional systematic uncertainty is 5.8% (see text for details). 



Decay Mode i 


Region j 


I 


II 


IV 


V 


VI 


VII 


1 i^*^(892)7r+ 


0.364 


0.360 






0.169 




2 scalar K*° 


0.352 


0.329 


0.357 


0.353 


0.372 


0.334 


4 p°{770)K+ 






0.287 


0.273 


0.286 




5 /o(980)K+ 








0.349 


0.362 




6 scalar / 


0.352 


0.329 


0.357 


0.353 


0.372 


0.334 


7 high mass 


0.352 


0.329 


0.357 


0.353 


0.372 


0.334 



Table 6: The branching fractions and uncertainties for the channels measured. 



Channel 


BF 

(xlO-*^) 


uncertainties (xlO ^) 


stat 


sys 


model 


interference 


K*^{892)tt+,K*^ K+TT- 


10.3 


±1.2 


±0.7 


+0.4 
-2.5 


±0.6 


"higher K*0"7r+,K*° ^ K+n- 


25.1 


±2.0 


±2.9 


+9.4 
-0.5 


±4.9 




184.6 


±3.2 


±9.7 






p°(770)K+,/)°(770) ^ vr+vr- 


3.9 


±1.2 


+0.3 
-0.6 


+0.3 
-3.2 


±1.2 


/o(980)K+,/o^7r+7r- 


9.2 


±1.2 


±0.6 


+ 1.2 

-1.9 


±1.6 


"higher fK+J ^ vr+vr- 


3.2 


±1.2 


±0.5 


+5.8 
-2.4 


±1.5 


Non-resonant 


5.2 


±1.9 


+0.8 
-1.8 


+3.3 
-7.5 


±6.4 




1.5 


±0.4 


±0.1 







4.2.1 Systematic Uncertainties 

In most regions, the largest contribution to the reconstruction systematic uncertainty, between 
3% and 11%, is from the uncertainties on the PDF parameters. The B background subtraction 
uncertainty is below 6% for all channels except for p^{770)K~^ (18%) and non-resonant 

(32%), where it is dominant. The systematic uncertainties from track efficiency corrections (2.4%), 
particle identification efficiency (4.5%), cosOt (2.8%) and the number of BB pairs (1.1%) are 
independent of the region considered. For "higher /" and "higher K*^" , there is an additional 
contribution of 9% to account for the selection efficiency being measured using scalar MC though 
the true angular distribution of the contribution is unknown. 

The model uncertainties include the effect of uncertainties in the resonance contributions, reso- 
nance masses, widths and line-shapes of the model, given in Table|31 The possibility that the dom- 
inant contribution to a region may be from another resonance is taken into account for Kq^{1A30) 
and /2(1270), where the alternative resonances are listed in Table |21 We consider the possibility 
that the component measured in Region VII does not extend into the other regions. We also allow 
for the possibility that the /o(980) resonance is described by the Flatte line-shape ^S] and the 
^0*^(1430) has the line-shape suggested by data from the LASS experiment ^2jUSl instead of the 
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Breit-Wigner form. We estimate model uncertainties on the branching fractions by varying the 
model contributions, resonance parameters and line-shapes within the parameters given in Table 01 
calculating a new P matrix and recalculating the branching fractions. The model uncertainty is 
the quadratic sum of all the variations in the branching fraction due to these individual changes. 

The effect of interference on the branching fractions is evaluated by generating many Dalitz 
plots with the observed branching fractions but with each contribution having a random phase 
and allowing interference, then measuring the branching fractions using the P matrix in the same 
way as done on the data. This produces a range of branching fractions for each channel with the 
correct branching fraction as the mean. The RMS variation is taken to be the uncertainty due to 
interference and this is shown in Table El under the column "interference" . 



4.3 Conclusion 

In conclusion, we have made preliminary measurements of the branching fractions for the fol- 
lowing channels with a statistical significance greater than 5o". This significance for a particular 
branching fraction Ba is evaluated by assuming that branching fraction is zero and minimizing, 
as a function of the other branching fractions, the separation of the measured yields and those 
obtainable from the branching fractions themselves. The significance is then given by \/x^{^a = 0). 

• B{B+ K*^{892)tt+,K*'^ K+tt-) = (10.3 ± 1.21^;?) x 10-^ 
. BiB+ ^ /o(980)K+, /o ^ TT+TT-) = (9.2 ± 1.2t|i) x IQ-^, 

• B{B+ XcoK+,XcO vr+7r~) = (1.46 ± 0.35 ± 0.12) x 10-^ 

• B{B+ D^i:+, D° K+TT-) = (184.6 ± 3.2 ± 9.7) x lO"*^ and 

• B{B+ ^"higher K*^''-k+) = (25.1 ± 2.0^^7°) x 10"^, where "higher means any combi- 
nation of K^°(1430),K2*°(1430) and i^i*°(1680). 

The first uncertainty is statistical and the second is systematic. This analysis has taken into account 
the uncertainty in the knowledge of the nature and parameterization of these resonances as well as 
interference between them, and these uncertainties are included (added in quadrature) in the above 
systematic value. Using the value of ^(K*° (892) -> K+tt") = 2/3, wefindi3(S+ ^ i^*0(892)7r+) = 
(15.5 lb 1.8]'^3'2) X 10~^, which is consistent with, and more precise than, previous measurements [TH 
I15| . The observation of the decay B'^ — > /o(980)if"'" is statistically significant, providing hints 
about the /o(980) production mechanism in the B system. There is also a significant signal for 
B~^ — > "higher i^**^" tt"*", which has a mass distribution which is partly described by the -ftrQ''(1430), 
as observed by the LASS experiment. 

We give 90% confidence-level upper limits for the branching fractions of the following channels, 
including the non-resonant component. This upper limit is taken as the value of that branching 
fraction for which the minimum separation of the measured yields and those calculated from 
the branching fractions is 1.64. 

• B{B+ p°(770)K+) < 6.2 X 10-^ 

• B{B^ A'+7r^7r+ non-resonant) < 17 x 10~^ and 

• i3(i?+ ^"higher f"K^) < 12 x 10~^, where "higher /" means a combination of /2(1270), 
/o(1370) and /2(1430). 
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The tight hmit on the non-resonant component means that it will be difficult to obtain the Unitarity 
Triangle angle 7 by the methods of Refs. ^ and P]. 
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Measurements of the Branching Fractions of Charged 
B Decays to iC+7r 7r+ Final States 

We present preliminary results of searches for exclusive charged- -B decays 
to K^TT^TT^ from 61.6 million BB pairs collected at the T{4S) resonance 
with the BaBar detector at the SLAC PEP-II asymmetric B Factory. The 
Dalitz plot is divided into eight regions and, using a maximum-likelihood 
fit, we measure statistically significant yields in all regions. We interpret the 
results as the following branching fractions averaged over charged-conjugate 
states: B{B+ K*^Ti+,K*^ K+tt-) = 10.3 ± 1.2^2? x lO^^, B{B+ 
foK+Jo ^ 7T+7T-) = 9.2±1.2l|i X 10-6, B{B+ ^ XcoK+,XcO ^ vt+tt-) = 
1.46 ± 0.35 ± 0.12 X 10-6 and B{B+ K+tt') = 184.6 ± 3.2 ± 

9.7 X 10-6. The first uncertainty is statistical and the second is system- 
atic and includes resonance-model and interference uncertainties. We give 
90% confidence-level upper limits on the branching fractions of the following 
channels: B{B~^ pK^) < 6.2x10-^ and B{B~^ K^-k~ti^ non — resonant) 
< 17 X 10-6. 



